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Abstract 


This  paper  extends  the  results  of  Penfield  and  Rubinstein  on  signal  delay  in 
RC  tree  networks.  Both  estimates  and  bounds  are  derived  for  the  step  response 
of  • leaky Y  RC  trees  and  meshes,  a  class  of  networks  that  is  appropriate  for 
modelling  interconnect  in  digital  bipolar  circuits.  This  paper  is  intended  to 
serve  as  a  tutorial  as  well  as  a  research  report.  Therefore  existing  work  is 
explained  in  some  detail  along  with  the  derivation  of  new  results. 
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I .  Introduction 


Simple  closed-form  bounds  for  signal  propagation  delay  in  linear 
RC  tree  models  for  interconnect  were  derived  in  [1].  The  theory  has 
since  been  extended  to  include  RC  mesh  models  as  well  [2].  The  goal  of 
this  report  is  to  further  extend  the  theory  to  allow  for  RC  networks  with 
grounded  resistors  (i.e.,  "leaky”  RC  trees  and  meshes  —  see  section  II 
for  a  precise  definition  of  these  terms) .  Both  delay  estimates  and  delay 
bounds  are  developed. 

The  practical  motivation  behind  allowing  resistive  paths  to  ground 
in  RC  mesh  circuits  is  an  attempt  to  more  accurately  model  interconnect 
on  bipolar  chips.  The  base  of  a  bipolar  tramsistor  loading  the  interconnect 
offers  a  (nonlinear)  resistive  path  to  ground.  A  corresponding  gate 
electrode  of  an  MOS  transistor  would  not  offer  such  a  resistive  path  to 
ground,  and  the  original  work  of  Penfield  et.  al.  relied  upon  the  absence 
of  such  paths. 


@  @ 
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Fig.  1:  Lumped  linear  RC  tree  model  for  MOS  logic  with  fanout 


Fig.  2:  Leaky  RC  tree  model  for  bipolar  logic  with  fanout 

Consider  the  RC  tree  in  Fig.  1.  The  voltage  source  e(t)  represents 
the  output  of  some  logic  gate.  The  output  of  this  gate  fans  out  through 
the  interconnect  to  the  inputs  of  two  other  logic  gates  at  nodes  @  and 
d)  •  Since  nodes  and  d)  represent  the  gate  terminals  of  MOS 
transistors,  there  is  no  need  for  a  resistive  path  to  ground  in  the  model. 

The  analogous  situation  in  bipoleu:  logic  is  illustrated  in  Fig.  2. 
Here  nodes  Q)  ©  represent  the  base  terminals  of  bipolar  transistors, 
so  we  need  a  resistive  path  to  ground  in  the  model.  The  resulting  circuit 
is  not  an  RC  mesh,  but  would  be  one  except  for  the  resistive  paths  to  ground 
This  is  the  type  of  network  discussed  in  this  report.  The  organization 
of  this  report  can  be  seen  in  the  table  of  contents. 


A  lumped  linear  ^  mesh  is  any  network  consisting  of  2-terminal 
linear  resistors  and  capacitors  driven  by  independent  sources ,  such 
that  one  side  of  every  capacitor  and  every  source  is  grounded,  no 
resistor  is  grounded ,  and  between  ^my  two  nodes  other  than  ground  there 
exists  a  path  consisting  entirely  of  resistors.  An  RC  mesh  is  called 
an  rc  tree  if  it  does  not  contain  a  closed  path  of  resistors. 

However,  in  this  paper  we  shall  consider  more  general  types  of 
circuits:  leaiky  RC  meshes  and  trees.  A  letOcy  RC  mesh  is  any  network 
meeting  all  of  the  requirements  of  an  RC  mesh  given  adsove,  except  for 
the  restriction  that  there  be  no  grounded  resistors.  So  for  every  leaky 
RC  mesh,  there  is  am  underlying  RC  mesh  which  is  obtained  by  open- 
circuiting  all  of  the  grounded  resistors  in  the  original  leaky  RC  mesh. 

A  network  is  a  leaky  RC  tree  if  it  is  a  leadcy  RC  mesh  whose  underlying 
RC  mesh  is  an  RC  tree  (i.e.,  there  are  no  resistor  loops  in  the  underlying 
RC  mesh) . 

III.  Network  Matrices 

In  this  section,  we  define  certain  network  matrices  to  be  used 
later  in  the  report.  Well-known  facts  about  these  matrices  will  be 
presented.  We  shall  motivate  these  definitions  by  the  example  circuit 
in  Fig.  3.  Note  that  this  example  circuit  would  be  an  RC  mesh  without 
the  grounded  resistor  R, . 


We  isolate  the  resistor  subnetwork  R  containing  all  the  resistors 


and  assign  reference  directions  to  the  capacitor  currents  i^,...,  i^^ 
as  shown  in  Fig.  4  {note  that  positive  current  is  defined  to  be  flowing 
out  of  the  capacitor  and  into  the  resistor  subnetwork) .  Let  the  node 
connected  to  the  independent  source  serve  as  the  datum  node  of  R.  Define 
the  matrix  R  as  follows : 

With  the  datum  grounded,  the  node  voltages  with  respect  to  ground 
are  given  in  terms  of  the  capacitor  currents  by  the  resistance  matrix  R 
as  shown  below: 


'  V^{t) 
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'  i^(t)’ 

i2(t) 

(1) 


We  can  introduce  vector  notation  to  write  (1)  in  a  more  compact  form. 

T 

Let  V  denote  the  column  vector  of  node  voltages  (v, ,  v. , . . . ,  v  )  .  Let 
.  1  2  N 

T 

1  denote  the  column  vector  of  capacitor  currents  (i.  ,  i^,)  ,  Then 

*  1  2  N 

(1)  becomes: 


v<t)  =  R  i(t) 


(2) 


So  with  the  datum  grounded,  the  numerical  value  of  any  matrix  entry  R^ ^ , 


in  ohms,  is  equal  to  the  potential  in  volts  we  would  observe  at  node  if 
a  1-amp  current  were  injected  into  node  ^  and  all  external  nodes  of  R 
other  than  ^  and  the  datum  were  open-circuited. 


In  the  special  case  of  an  RC  tree  the  matrix  entries  can  be  read 
off  by  inspection,  since  in  that  case  is  simply  the  sum  of  the  resistances 
along  the  route  obtained  by  intersecting  the  path  from  node  ^  to  datum 
with  the  path  from  node  to  datum.  So,  the  topological  definition 
of  R^^  for  an  RC  tree  given  in  [1]  is  consistent  with  the  more  general 
interpretation  of  R^^  as  an  element  of  the  resistance  matrix  R  describing 
the  resistor  subnetwork  R. 

From  circuit  theory,  the  resistance  matrix  R  has  certain  well-known 
properties  independent  of  the  particular  resistor  subnetwork  being 
described  by  R: 

(i)  R  is  symmetric  (i-e.,  R^^  =  R^^  for  all  i,j):  This  is  a 

consequence  of  R  being  a  "reciprocal"  network.  [3;  Chapter  16] 

(ii)  R  is  positive-definite  (assuming  all  the  resistors  in  R  are 

positive) :  This  is  true  because  the  power  dissipated  in  R, 

T 

given  by  the  quadratic  form  i  R  i  is  positive  for  any  nonzero 
choice  of  currents  i. 

These  properties  of  R  are  quite  standard  and  play  only  a  minor  role  in 
t.he  theory.  But  R  possesses  more  structure  than  symmetry  and  positive¬ 
definiteness,  and  this  additional  structure,  described  in  Section  VI  and 
proved  in  Appendix  A,  is  the  foundation  for  the  theory  developed  in  this 


report. 


b.)  The  conductance  matrix  G: 


data’Tt 


Fig.  5:  The  resistor  subnetwork  K  extracted  from  the  circuit 
in  Fig.  3  (driven  with  voltage  sources) 


Again,  we  isolate  the  resistor  subnetwork  R  (Fig.  5) .  With  the 
datum  grounded,  we  now  express  the  capacitor  currents  in  terms  of  the 
node  voltages  by  the  conductance  matrix  G  as  shown  below: 


ij^  (t) 

i2  (t) 


i^(t) 


Gii  G^2  •••  ^iN 

^21  ^22  '^2N 


G  G  .  .  .  G 
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Vj^(t) 
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or,  using  the  vector  shorthand: 


i  (t)  *  G  v(t ) 


So  with  the  datum  grounded,  the  numerical  value  of  einy  matrix  entry 
in  m.nos  ,  is  equal  to  the  current  in  amps  we  would  observe  flowing  into 


3 


node  (T)  if  a  1-volt  source  were  applied  to  node  and  all  external  nodes 
of  R  other  than  were  grounded. 

If  R  is  the  resistance  matrix  for  the  same  resistor  subnetwork 
described  by  G,  then  G  =  R  Like  R»  G  is  also  positive-definite 
symmetric  (for  similar  reasons  as  R) /  and  G  possesses  additional  structure 
described  in  Appendix  A  which  is  crucial  to  the  results  in  Section  VI. 

c . )  The  capacitance  matrix  C : 

This  matrix  is  not  as  central  to  the  theory  as  either  the  R  or 
the  G  matrix.  It  is  simply  notational  shorthand  for  writing  down  the 
capacitors  in  the  RC  network  under  consideration.  We  define 
I  C  =  diagCCj^,  C^}.  So  C  is  a  “diagonal"  matrix,  (i.e.,  has  off- 

diagonal  elements  equal  to  zero) ,  eUid  the  diagonal  elements  are  equal 
to  the  capacitor  values  consistent  with  the  assigned  node  numbers. 

1  One  simple  property  to  note  about  C  is  that  its  inverse,  C  is 

also  diagonal:  C  ^  =  diagfc^  (Of  course,  C  ^ 

doesn't  exist  in  the  degenerate  case  where  some  individual  capacitors 
I  are  equal  to  zero.) 

r  IV.  Network  Differential  Equations  in  Ccuionical  Form 


In  order  to  write  down  the  network  differential  equations,  we  must 
first  coneider  what  happens  when  the  datum,  e(t),  is  not  grounded.  In 
the  previous  section  we  wrote  the  node  voltages  in  terms  of  the  capacitor 
currents  when  the  datum  was  grounded: 


v(t)  *  R  i(t)  ;  (dat'jn  e(t)  =  0)  (5). 

When  e(t)  is  not  zero,  then  by  superposition  we  can  write: 

v(t)  =  R  i (t)  +  xe (t)  (6) , 

T 

where  x  =  (x, ,  X-,...,  x^,)  is  a  column  vector  of  dimensionless  numbers 
^12  N 

(each  x^  satisfying  the  inequality:  0  <  x^  £  1) .  Referring  back  to  Fig.  4, 
each  is  numerically  equal  to  the  potential  (in  volts)  produced  at  node 
^  due  solely  to  a  1-volt  source  at  the  datum,  while  open-circuiting  all 
of  the  external  current  sources. 


Example In  Fig.  3,  if  R^  =  ®  (open-circuit),  the  circuit  is  an  RC 

T 

mesh  and  x  =  (1111)  . 

Example:  In  Fig.  3,  if  R^  is  some  finite  resistance,  x  requires  some 

^  ■  b  ^ 

computation.  The  result  is: 


(R4+R5) 


R-  +  (R.+R,)  II  (R.+R.) 


R,(R.+R,) 

R  +  -.3  5 

6  R,+R,+R  +R, 


2  3  4  5 


'6  R-+R,+R  +R^ 

2  3  4  5 


The  general  computational  problem  of  calculating  x  is  discussed  in  Section 


VIII. 


To  begin  the  formulation  of  the  network  differential  equations, 
we  start  by  rewriting  equation  (6), 

X  e(t)  -  v{t)  =  -R  i(t)  (7) . 

The  constitutive  relations  for  the  capacitors  can  be  neatly  written  as: 

i (t)  »  -C  v(t)  (8) , 

where  the  minus  sign  arises  from  our  sign  convention  for  capacitor  currents. 
Substituting  (3)  into  (7)  we  have, 

X  e(t)  -  v(t)  =  [R  C]  v(t)  (9). 

The  system  of  equations  in  (9)  is  the  canonical  form  for  leaky  RC 
meshes.  The  equations  in  (9)  are  also  the  most  general:  valid  for  all 
time  and  for  any  possible  e(t).  At  this  point,  it  is  worthwhile  specializing 
the  system  in  (9)  somewhat  by  considering  only  a  certain  form  for  e(t). 

Since  the  voltage  e(t)  models  the  output  of  a  logic  gate  undergoing 
a  logic  transition,  typically  e(t)  reaches  eund  remains  at  its  final  value, 
e(«>),  for  some  time  interval  t>T. 


At  equilibrium  v(t)  =  0,  so  (10)  becomes: 


xe(»)  -  V  =  [RC]0  =  0 
.eq  --  - 


V  =  xe(®) 
-eq 


(11) 


Finally,  since  estimates  for  and  bounds  on  the  step  response  will  be 
derived  later  in  this  report,  the  Ccmonical  equations  for  a  step  input 
will  be  given.  Using  (10)  and  (11) , 


Network  Differential  Equations  for  a  Step  Input: 


V  -  v(t)  “  [RC]v(t)  ;  t>0 
,eq  -  - 


(12) 


V.  Delay  Estimation 


In  this  section,  a  good  first-order  estimate  of  the  delay  through 


a  leaky  RC  mesh  is  developed.  This  delay  estimate  is  called  the  "Elmore 


Delay"  in  the  literature  [4],  also  [1],  [5-7].  We  will  denote  it  here  by 

"T^  ".  By  convention,  we  call  the  output  node  of  interest  "node 
i 

So  the  subscript  on  simply  indicates  that  the  delay  estimate  depends 

i 

on  which  node  in  the  RC  network  we  have  chosen  as  our  output  node.  Finally, 


since  is  just  a  single  number,  the  use  of  alone  takes  no  account 

of  the  threshold  voltage  of  the  device  driven  by  node  .  A  simple 


procedure  for  accounting  for  the  logic  threshold  is  presented. 


a.)  Definition  of  T_  • 

- D. 

1 


We  define  T  in  terms  of  the  response  at  the  output  node  (T  to 

r\ 


a  step  transition  at  the  input.  We  allow  for  a  non-zero  initial 


I 


equilibrium  at  the  time  of  the  step  input,  and  denote  the  initial  equilibri'iir. 

T 

voltage  distribution  with  the  vector  v^  =  (v^  ,  v^  ,  - . .  ,  v  )  .  Using 

o  o  o 

the  same  notation  as  the  previous  section,  we  denote  the  final  equilibrium 


voltage  distribution  with  the  vector  v 
- ^ - .eq 


(v  ,  V  , . . . ,  v„  )  .  The 
1  ^  N 

eq  eq  eq 


step  transition  can  be  either  a  step  increase  (v  >  v  )  or  a  step  decrease 

,eq  ,o 

(v  <  V  ) .  The  amplitude  of  the  step  input  is  equal  to  —  (v  -  v  ) . 

-eq  -o  — ^ ^  X.  1  1 

1  eq  o 


Definition: 


“i  -  ''1'’=' 
eg 


V.  -  V. 

1  1 

eq  o 


where  v^ (t)  is  the  response  at  node  to  a  step  transition  at  the  input 


at  time  t=0  of  amplitude  —  (v.  -  v.  )  . 

X .  1  1 

1  eq  o 

Graphically,  T^^  is  shown  below: 


''i  - 

eq 


V  -  V, 

eq  o 


shaded  area  =  T 


Definition  of  T 


b. )  Computation  of  T_  : 


Rewriting  the  canonical  matrix  equation  for  a  step  input,  we  have 


V  -  v(t)  =  (RC]v(t)  ;  t>0 
-eq  -  --  - 


Taking  the  i  row  of  this  matrix  equation, 


eq  k 

(where  J  is  taken  over  all  capacitor  nodes  in  the  RC  network) . 
k 

Dividing  by  v.  -  v.  and  then  integrating  from  t=0  to  t=*,  we  obtain 
^eq  o 

•  OO 

D.  ■  “  _  r 


.  V.  -  V.  (t) 

All 


V.  -  V. 

1  1 
eq  o 


.  IX  K  ; 


V  (t)  dt 

k 


V.  -  V. 

1  1 
eq  o 


I  ■  \  > 

k  eq  o 


V,  -  V. 

1  1 
eq  o 


So  in  order  to  compute  we  need  to  know: 

i 

1)  The  initial  and  final  equilibrium  voltage  distributions,  v^ 

and  Note  that  this  is  really  equivalent  to  knowing  the  single 

vector  X,  since  for  a  step  input  we  have  v^  =  xe(0  )  and  v^^  =  xe (0  ). 

In  fact,  we  can  express  T^  in  terms  of  x  instead  of  v  and  v 

D .  ,  -o  -eq 

1 


T„  =  —  y  R-i-C,  X 
D.  X.  “  ik  k  k 
1 


2)  The  1  row  (or  column)  of  the  resistance  matrix  R. 

3)  The  capacitance  matrix  C  (this  requires  no  computation) . 

Note  that  when  we  are  dealing  with  an  RC  mesh  (i,e.,  no  resistive 

paths  to  ground)  starting  from  a  zero  initial  equilibrium,  then  v  =0  and 

~  o  ' 

V,  =  V,  for  all  k7<i.  Hence  T_  reduces  to  7  R.,  C,  ,  which  is  consistent 
1  k  D.  V  k 

eq  eq  i  k 

with  the  definition  given  in  earlier  works  [1] ,  [2] . 


c.)  Alternate  Interpretation  of  T, 


We  have  already  defined  T^  in  terms  of  the  response  observed  at 

node  ^  ,  v^(t),  to  a  step  input  applied  at  time  t=0,  when  the  RC  network 

is  at  some  (possibly  non-zero)  equilibrium.  Alternately,  we  can  interpret 

T  in  terms  of  v. (t) .  Note  that  v. (t) ,  in  addition  to  being  the  time 
D.  1  1 

1 

derivative  of  v^(t)  is  also  the  response  observed  at  node  ^  to  an  impuls 
input  applied  at  time  t*0,  when  the  RC  network  is  at  zero  equilibrium. 

We  will  show  in  the  next  section  that  the  step  response  v^Ct)  is 
monotone,  which  means  that  v^(t)  can  be  interpreted  as  a  density  function 
for  a  continuous  random  variable,  when  properly  normalized: 


(y, (t) )  1  •  j 

normalized 


V.  (t) 
1 


r- 

1 

.  o 


.  (t)dt 


v^(t) 

V,  -  V. 

1  1 
eq  o 


T_  is  the  first  moment  (or  "center  of  mass")  of  this  normalized  zero- 


state  impulse  response,  as  shown  below 


Integrating  by  parts, 


where  we  have  used  the  fact  that  limit  t(v.  -  v. (t))=0  because  v, (t)  -►  v. 

t--»  eq  eq 

exponentially  as  t-«“. 


d. )  Incorporating  Logic  Threshold: 


As  mentioned  earlier,  the  use  of  alone  takes  no  account  of  the 

i 

threshold  voltage  of  the  device  driven  by  node  Q  .  To  do  so  we  need  an 
estimate  for  the  entire  step  response  voltage  waveform  v^(t).  A  reasonadsle 
estimate  is 

V.  (t)  =  V.  +  (v.  -  v^  )  e  ^  (20) 

eq  o  eq 

since,  like  the  exact  but  unknown  response: 


1)  V. (t)  starts  at  V.  . 

2)  v^(t)  is  monotone. 

3)  V. (t)  asymptotes  to  v. 

"*  Vi  -  v^(t) 

v-^-  V—  = 
•'Oil  1 

eq  o 


In  the  next  section  we  will  show  that  v^{t)  always  lies  between  the  best 
currently  available  upper  and  lower  bounds  on  the  exact  step  response  v^(t) 
We  can  invert  this  estimate  of  the  entire  waveform  in  order  to 
estimate  the  time  to  reach  a  given  logic  threshold  voltage  v  : 


V.  (V,  -V.  )e 

1  11 
eq  o  eq 


t  =  In 
D. 

1 


V.  -  V. 

1  1 

eq  o 

V-  -  V 
1  T 

eq 


(21)  . 


Example:  If  the  logic  threshold  of  the  device  driven  by  node  (T)  is 

v_  »  V  1  (V.  -  V.  ),  then  instead  of  using  T_  as  a  delay  estimate 


we  use  the  somewhat  lower  estimate,  in  2  *  0.69  . 

i 

Example:  If  the  logic  threshold  of  the  device  driven  by  node  is 

V  =  V  +  (1  -  -)  (V.  -  V.  )  =  V.  +  0.63  (v.  -  V.  ),  then  we  use 

Ti  eiii  11 

o  eq  o  o  eq  o 

exactly  as  a  delay  estimate. 

i 

The  estimate  (20)  is  not  generally  exact:  for  one  thing,  it  contains 


only  a  single  natural  frequency:  - 


To  better  appreciate  the  nature 


of  the  approximation  involved,  Horowitz  [7]  has  noted  that,  in  the  case 
of  no  grounded  resistors ,  a  sufficient  condition  for  (20)  to  be  exact  is 
that  the  time  derivatives  of  all  the  capacitor  voltages  are  equal:  i.e., 
Vj^(t)  *  v^(t)  for  all  }c.  In  the  case  where  we  ^  have  grounded  resistors , 
the  estimate  (20)  is  exact  when  the  normalized  time  derivatives  of  all  the 
capacitor  voltages  are  equal: 


v^(t) 


V.  (t) 
1 


V,  -  V,  V.  -  V. 

)(  k  1  1 

eq  o  eq  o 


for  all  k?<i. 


(22)  . 


We  can  see  this  by  going  back  to  the  i  row  of  the  original  system  of 
equations  for  the  node  voltages: 


V.  -  V. (t) 
1  1 
eq 


I  'ik'k 


k  k 

^  ^  °  v^(t)  ;  using  the  approxima* ion 

in  (22) 


I  ■  \  ' 

eq  o 


V.  -  V. 

1  1 
eq  o 


V.  (t) 
1 


T_  V. (t) 


(23)  . 


The  solution  to  the  first-order  differential  equation  given  by  (23)  is 


v^  (t)  =  (t) . 

Another  view  of  the  nature  of  the  approximation  (20)  will  emerge 
from  consideration  of  the  step  response  bounds,  which  we  consider  next. 

VI .  Optimal  Control  Method  for  Determing  Bounds  on 
the  Step  Response 

The  purpose  of  this  section  is  to  clearly  derive  step  response  bounds 
for  leaky  RC  meshes  using  an  optimal  control  point-of-view.  The  bounds 
derived  will  be  analogous  to  those  in  [1].  Recent  bound-tightening 
results  such  as  slew-rate  limitations  [8] ,  and  (in  the  case  of  RC  trees) 
spatial  convexity  [9]  will  not  be  considered  here.  These  recent  improvements 
are  omitted  because  they  do  not  add  to  the  basic  understanding  of  the 
optimal  control  approach.  We  do  not  yet  know  exactly  to  what  extent 
these  bound- tightening  results  cam  be  applied  to  leaky  RC  meshes.  It  is 
clear,  however,  that  at  least  the  simplest  of  the  slew-rate  limitations 
applies  and  always  gives  somewhat  tighter  bounds  than  those  developed  here. 

a. )  Preliminary  Results; 

Fact  1;  For  any  three  nodes  0,  (2)  »  0  of  a  leaky  RC  mesh, 

R.  ,R.  .  >  R.  .R.  . 

11  K3  —  ki  13 

where  the  resistances  in  (24)  aure  elements  of  the 
resistance  matrix  R  in  (1). 


The  proof  of  Fact  1  depends  on  certain  properties  of  the  conductance  matrix 
G.  These  properties,  as  well  as  the  proof  of  Fact  1,  are  given  in 


Appendix  A. 


Fact  2 : 

For  any  node  ^  of  a  leaky  RC  mesh,  the  observed 

response,  v^ (t) ,  to  a  step  transition  at  the 

input  applied  at  time  t=0,  when  the  network  is 

at  some  (possibly  non-zero)  equilibrium,  is 

monotone : 

V. (t)  >  0,  Vt  >  0  ;  if  V.  >  V, 

^  ^eq 

V. (t)  <  0,  Vt  >  0  ;  if  V.  <  V. 

^  ~  ■  'eg  ^o 

(25)  . 

The  proof  of  Fact  2  is  also  given  in  Appendix  A. 


Fact  3:  For  any  two  nodes  of  a  leaky  RC  mesh: 

On  an  "up"  transition  (v  >  v  )  we  have, 

.eq  .o 

eq  eq 

eq  eq 

On  a  "down"  step  transition  (v  <  v  )  the  sense  of 

.eg  -O 

each  inequality  is  reversed. 


eq  eq 


The  proof  of  Fact  3,  which  uses  both  Fact  1  and  Fact  2,  is  presented  here 


in  the  body  of  the  report. 


Proof  of  Fact  3 : 


The  system  of  differential  equations  for  the  node  voltages  is: 


V  -  v(t)  =  [RC]v(t) 
,eq  -  --  - 


till  P 

Taking  the  k  row  -  v(t)  =  I  R  .C.v.(t) 

K  K  .  ^  J  T  T 

eq  3 

til  r  . 

Taking  the  i  row  =>  v.  -  v.(t)  -)  R..C,v.(t) 

eq  ,  j  j  j 


(30) 


(31) 


(32) 
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Using  (31)  and  (32)  ,  we  compute: 


-  ''k""  -  '>ki'\  -  Vjit  ))  -  I  (R 

eq  eq  j  j  j  j  j 


(33) 


By  applying  Facts  1  and  2,  (26)  and  (28)  now  follow  directly  from  (33). 
By  interchanging  the  subscripts  i  k  and  using  the  fact  that  R  is 
symmetric,  (27)  and  (29)  follow  from  (26)  and  (28)  respectively. 


b.)  A  New  State  Variable: 


At  this  point,  we  will  introduce  a  new  state  varicible,  f^,  that 
will  be  part  of  the  reduced  order  model  in  the  formulation  of  the  optimal 
control  problem.  We  also  will  prove  a  lemma  concerning  f^.  The  result 
of  this  lemma  will  serve  as  the  state  constraints  for  the  optimal  control 
problem.  The  presentation  of  the  optimal  control  problem  itself  is 
delayed  until  later  (part  (c)  of  Section  VI). 

The  new  state  variable  is  defined  as  follows: 


f^(t) 


“>  V,  -  V.  (f  ) 
1  1 

eg _ 

V  ,  -  V 

^  1  1 

t  eq  o 


(34)  , 


where  (as  usual):  f  v^(t)  is  the  response  at  node  to  a  step  input 


at  time  t=0. 


V.  is  the  initial  equilibrium  voltage  at  node  Q) 
o 

V.  is  the  final  equilibrium  voltage  at  node  . 
^eq 


This  definition  is  illustrated  below: 


Note  the  similarity  of  Fig.  7  to  Fig.  6  (where  we  defined  .  In  fact, 

i 

we  can  verify  directly  from  these  illustrations  the  initial  condition: 


f .  (0)  =  (35)  . 

1  D. 

1 


Due  to  normalization  in  the  denominator  of  (34),  the  behavior  of  f^(t)  is 


the  same  for  "up"  steps  and  "down"  steps  of  arbitrary  magnitude  at  the 
input.  Starting  at  T  when  t=0,  f . (t)  will  decrease  monotonically  towards 


Proof  of  Lemma: 


Starting  with  the  i  row  of  the  canonical  matrix  differential 


equation, 


eq  k 


(39)  , 


we  divide  by  v,  -  v.  ,  and  then  integrate  from  t’=t  to  t'=®  to  obtain 
^eq  ^o 


f ,  (t)  : 
1 


f .  (t,  ^  -S3 -  df  =  Is - 

1  V .  -  V  ,  V ,  -  V  , 

^11  11 
•'t  eq  o  eq  o 

^  R.,  C,  (V,  -  V,  (t)  ) 

,  ik  k  )c  k 

k  eq 


V.  -  V , 

1  1 
eq  o 


(40)  . 


Using  (26)  from  Fact  3  on  an  "up"  step  (or  (28)  on  a  "down"  step]  in  (40) 


gives : 


)  ik  k 


^  (V.  -  V, (t) 


where  we  have  used  the  fact  that  by  symmetry  of  R. 

Using  (27)  from  Fact  3  on  an  "up"  step  [or  (29)  on  a  "down"  step] 


in  (40)  gives: 


I  R-vC, 


f.  (t)  < 
1  — 


r  'Sck  1 

*-  ki  eq  -> 


■ 

k 


V.  -  V. 

1  1 
eq  o 

/V.  -  V.  (t) 

1  1 

—§3 _ 

V.  -  V, 

I  1 

'■  eq  o 


(42)  , 


where,  again,  we  have  used  the  fact  that  R^^^  =  R^^^.  Together,  (41)  and 


(42)  constitute  (36) .  This  proves  the  lemma. 


c.)  The  Optimal  Control  Approach: 


The  step  response  bounds,  as  in  [1],  can  be  derived  by  "ad  hoc" 
manipulation  of  (25)  and  (36),  but  the  methodology  is  somewhat  obscure. 

A  clearer  view  emerges  from  recasting  these  calculations  into  the  form  of 
a  linear  minimum  -  (and  maximum  -)  time  optimal  control  problem  with  state 
constraints.  In  addition,  this  approach  has  the  advantage  that  it  becomes 
much  clearer  how  to  incorporate  additional  information  in  order  to  obtain 
tighter  bounds. 

The  optimal  control  approach  is  based  on  a  reduced  order  model  of 
the  voltage  transition  process  at  a  selected  output  node  .  We  regard 
this  process  as  being  described  by  the  following  second  order  model: 


where : 


(i)  g^(t)  is  shorthand  for  the  normalized  "distance  to  go"  for 


the  step  response  at  node  ^  : 


V,  -  v.(t) 
i.e.,  g. (t)  = — ^3 - 

1  V.  -  V, 

1  1 
eq  o 


(ii)  f^(t)  is  defined  as  in  (34). 

^  V.  -  V.  (f  ) 

A  r“  1  1  A  f* 

[Note  that  f^(t)  =  ^  ^ -  dt'  =  g.(t')dt’, 

J  t  ^  i  J  ^ 

so  by  the  fundamental  thISrem  of  calculus  we  have  f^(t)  = 

-g^ (t) ,  as  in,  (43) .  ] 

(iii)  The  input  u(t)  is  introduced  to  represent  the  unknown  waveform 
r  V,  (t)  ^ 


g.  (t) 


V,  -  V  . 

1  1 
eq  o 


In  addition  to  the  dynamics  in  (43)  we  also  have, 
initial  conditions  (see  (35)):  F  F 


9^(0) 


state  constraints  (see  (36)): 


g^(t)  <  f.  (t)  <  Tpg^(t) 
i 


and  the  input  constraint  (see  (25)): 


(45)  , 


u(t)  <  0 


To  arrive  at  upper  and  lower  bounds  on  g^(t)  we  pose  the  following 

* 

optimal  control  problem:  What  is  the  minimum  time  t  .  (g.)  and  the  maximum 

min  1 

* 

time  t  (g.)  required  for  a  solution  of  (43)  to  pass  through  a  given  "target 
insix  ^ 

It  * 

value  g^  (with  0  ^  1)  amd  still  satisfy  the  constraints  given  by  (45) , 

*  * 

(46),  and  (47)?  Once  the  maximum  (minimum)  time  t  (g.)  (t  .  (a.))  is 

max  1  min  ’i 

* 

Icnown  as  a  function  of  g^,  then  the  inverse  function  gives  am  upper  (lower) 
bound  on  g^(t).  [An  upper  (lower)  bound  on  g^(t)  can  then  be  easily 
converted  to  a  lower  (upper)  bound  on  the  voltage  step  response  v^(t).] 

The  above  problem  is  called  a  "state  constrained  linear  optimal 
control  problem"  and  the  Pontryagin's  maximum  principle  could  be  used  to 
find  the  optimal  u(t)  amd  the  resulting  optimal  state-space  trajectory 
(f^(t),  g^(t)).  However,  the  problem  cam  be  solved  in  a  much  simpler  way 
as  follows.  Regard  amy  trajectory  in  the  plame  as  a  function  of 

f^,  i.e.,  g^(f^).  From  (43)  we  have 


dt  =  -  —  df . 

9:  i 


(48)  , 


amd  the  time  required  by  the  trajectory  to  pass  through  two  points 


(f .  (0) ,  g.  (0))  and  (f.(T),  g.(T))  is 


(f^(0) ,g^(0) ) 


slope 


Fig.  8:  "Typical”  trajectory  in  the  plane 


The  trajectory  corresponding  to  a  voltage  transition  at  node  ^ 
starts  at  (f.  (0)  ,  g^(0))=  (TJ^ ,  1)  and  because  f .  (t)  _<  0  and  g.  (t)  ^  0,  it 
moves  down  emd  to  the  left  towards  the  origin  (see  Fig.  8) .  The  trajectory 
is  confined  to  the  region  between  g^  »  '^i  *  T^^i  stated  by 

the  state  constraints  (46).  From  (49),  it  can  be  seen  that  for  g^(t) 
to  change  from  1  to  a  particular  value  g^,  the  lowest  (highest)  possible 
trajectory  (i.e.,  with  the  minimum  (maximum)  possible  g^(t))  over  the 
widest  (narrowest)  integration  interval  on  the  f^-axis  requires  maximum 
(minimum)  time  and  therefore  is  the  slowest  (fastest)  trajectory.  We  can 


then  calculate  t  (g.)  and  t  .  (g.)  for  the  optimal  trajectories,  using 
max  min  i 


Solution  to  the  maximum- time  problem: 


There  are  two  different  algebraic  forms  for  t  (g.)  over  the  range 

mdx  1 


0  <  g.  <1.  This  is  a  result  of  the  fact  that  the  maximum-time  state  _ 

1  -  T„ 

*  Dj_ 

space  trajectory  assumes  a  more  complex  form  as  g.  is  decreased  below  - — 

^  T  ^ 

D.  * 

Case  #1  ^  1  1  ^  ■ 

P 


using  (49)  , 


df . 
1 


t 

max 


(gi) 


^  df .  + 

T„  *  g. 

R.g.  1 
11 


T„g- 

Pi 


* 


* 

R.g. 

11 


df.  +  T 
1  P 


T 

r  D, 


T„g- 

Pi 


*  [Tpg^  -  +  Tp  [in(Tp  )  -an(Tpg*)l 


'vl 

T  -  T_  +  T„in 

P  R.  P 

1 

* 

T„g. 

1  Pi; 

Solution  to  the  minimum-time  problem; 

* 

There  are  three  different  algebraic  forms  for  t  .  (g.)  over  the 

min  1 

range  0  g^  £  1*  This  is  a  result  of  the  fact  that  the  minimum-time 

state  space  trajectory  assumes  a  progressively  more  complex  form  as 

•  1  i 

g.  is  decreased  below  - —  and  - —  . 


Fig.  11:  Minimum-Time  Trajectory  (for  ^  9-  ^  D 

P 

Using  (49) ,  we  see  that  there  is  no  area  (along  the  f^-axis)  under 
the  minimum- time  trajectory,  so: 


The  results  in  (50)  -  (54)  can  be  suznmarized  in  the  following  graph  of 

#  *  *  * 

both  t  (q  )  and  t  .  (g.)  vs.  g.  (over  the  entire  range  0  <  g.  <  1) : 
max  i  min  i  i  i  — 


* 


Fig.  14:  Maximum  and  Minimum  Times  as  a  Function  of  g^ 


Since  g^(t)  ^0,  at  any  instant  of  time  t  the  actual  value  of 

* 

g^(t)  is  less  than  or  equal  to  the  value  of  g^  that  is  reached  via  the 
maximum-time  trajectory  in  time  t,  and  is  greater  than  or  equal  to 
the  value  of  g^  that  is  reached  via  the  minimum-time  trajectory  in 

It 

time  t.  For  this  reason,  the  inverse  functions  of  t  ,  (g. )  and  t  .  (c 


Finally,  the  bounds  on  g^(t)  can  be  easily  converted  to  bounds  on  the 

voltage  step  response  v^(t)  as  follows: 

'  v^  -  V^(t)' 

2.(t)  <  g.(t)  =7-^2— -  <g.(t)  (58) 

i  i 
eq  o  ■' 


V.  (t)  <  V.  (t)  <  V.  (t) 

—  1  ~  1 

where , 

* 

V.  t)  =  V.  -  (v.  -  V.  )g. (t) 

—1  1  1  11 
eq  eq  o 

V. (t)  “V.  -  (v,  -  V.  )  g. (t) 

1  1  1  1  -^1 
eq  eq  o 

These  bounds  are  plotted  below: 


For  any  leaiky  RC  mesh  it  is  always  true  that  <  T  <  T  (evaluate 

R.  —  D.  —  P 
1  1 

the  state  constraints  (46)  at  time  t=0;  see  also  Appendix  B) .  The  Elmore 

waveform  estimate,  v^(t),  represents  em  effort  to  approximate  the  dynamics 

of  a  higher  order  network  by  one  with  a  single  time  constant  T  :  the 

»  i 

estimate  is  exact  only  in  that  case.  However,  v^(t)  is  a  very  good 
estimate  whenever  T  -  T  T  because  this  means  the  state-constraint 

&  K .  u . 

1  1 

cone  shown  in  Fig.  8-13  is  very  n^u:row,  so  the  actual  trajectory  in  the 
f^-g^  plane  is  well-approximated  by  the  Elmore  trajectory  of  a  straight 
line  of  slope  .  [In  fact,  the  voltage  bounds,  v^(t)  and  v^(t),  are 
very  tight  whenever  T  -  T  <<  T  for  exactly  the  same  reason.  The 

<r  I\  ,  D  . 

1  1 

minimum-  and  maximum-time  trajectories  are  very  close  together  simply  by 
virtue  of  the  fact  that  they  each  must  lie  inside  a  very  narrow  state- 
constraint  cone.] 


VII.  Bounds  for  Non-Step  Inputs 

In  the  previous  section  we  derived  rigorous  upper  and  lower 
bounds  on  the  step  response  voltage  waveform  at  any  selected  output  node 
A  question  which  naturally  arises  is:  How  do  we  bound  the  response 
at  node  to  an  arbitrary  non-step  input?  In  this  section,  we  will  show 
a  general  method  for  deriving  bounds  for  an  arbitrary  input  in  terms  of 
the  step  response  bounds.  We  will  first  consider  the  case  where  the 
leaky  RC  mesh  is  at  zero  initial  equilibrium  (this  treatment  is  essentially 
the  Seune  as  the  one  given  in  Appendix  E  of  reference  [1]).  We  will  then 
show  how  these  results  can  be  easily  adapted  to  an  arbitrary  (possibly 
non-zero)  initial  equilibrium.  Finally,  we  will  carry  out  the  calculations 


Li.  u  j*  jj 
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explicitly  for  a  particular  example:  a  saturating  ramp  input, 
a. )  Non-Step  Input  (Zero  Initial  Equilibrium) : 

The  method  presented  here  depends  on  the  fact  that  the  RC  network 

is  linear  (and  time-invariant).  We  will  first  introduce  some  notation: 

A 

h^(t)  =  zero-state  response  at  node  to  a  unit  impulse  input. 
w^(t>  *  zero-state  response  at  node  ^to  a  \init  step  input;  so 
w^ (t)  =  h^ (t) . 

[Here  we  use  the  notation  w^(t)  rather  than  v^(t).  The  reason 

for  this  is  that  w^(t)  denotes  a  more  specific  waveform  than  was  denoted 

by  our  previous  usage  of  v^(t).  We  had  previously  used  v^(t)  to  represent 

the  response  to  a  step  input  of  arbitrary  magnitude  starting  from  some 

arbitrary  initial  equilibrium.  Note  also  the  difference  in  dimensions: 

v^(t)  has  dimensions  of  volts,  but  w^(t)  is  dimensionless  since  it 

represents  the  unit  step  response.) 

A 

e(t)  =  arbitrary  non-step  input. 

y^(t)  =  zero-state  response  at  node  to  the  input  e(t). 

Using  Laplace  treinsforms, 

W,  (s)  =  H.  (s)  ■  -  (61)  . 

1  IS 

Y^(s)  =  H^(s)  ■  E(s)  (62) . 

. 

Combining  (61)  and  (62)  we  have, 

Y.  (S)  *  [-^.(s))  [sE(s)l  -  (W  (s)I  [sE(s)]  (63). 

1  s  1  1 

But,  sE(s)  »  ^de (t)  ,  _  assuming  e(0)=0 


(64)  . 
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So  converting  (63)  back  to  the  time  domain  (and  using  (64)),  we  have: 


w. (t) * 
1 


de(t) 

dt 


■t 

w  (t-f) 

0 


de(t') 

dt 


df 


(65)  . 


The  step  response  bounds  for  w^(t)  are  given  by: 


w^(t)  ^  w^(t)  w^(t) 


(66)  , 


where  (referring  back  to  (59))  we  have: 


w^(t)  =  x^(l-g^(t) ) 


w^(t)  =  X^(l-£^(t)) 


(67)  , 


since  w.  =0  and  w,  *  x. . 
1  11 
o  eq 


In  order  to  obtain  bounds  on  y^(t),  we  note  that  convolution  of 


(66)  with  retains  the  sense  of  the  inequality  when  >  0  Vt>0, 

at  Qw  "" 


and  reverses  the  sense  of  the  inequality  when  <  0  Vt>0,  i.e., 


dt 


de(t  ,  ,  de  t)  >  ,  -  de  t  de(t)  ^ 

w,  (t  *  — r- —  <  (y.  (t)  =  w.  (t)  *  — TT —  )  <  w,  (t)  *  —  ;  —rr —  >0  Vt  > 

—1  dt  —  1  1  dt  —  1  dt  dt  —  — 


w.  (t)  *  <  (y.  (t)  -  w,  (t)  * 


de(t) 


dt  -  -"i 


dt 


)  1  w.  (t)  *  <0  Vt  >  0 

dt  dt  ~  — 


- ^1 


(68)  . 


So,  when  the  arbitrary  input  e(t)  is  monotone ,  the  bounds  on  y^(t)  are 


given  directly  by  (68) : 


Zi(t)  = 


y^ft) 


. de ( t ' )  , .  .  . 

w^(t-t')  —  ;  e(t)  increasing 


-  de(t') 

w.{t-t')  . -, —  ;  e(t)  decreasing 

J  0  ^ 

w.(t-t')  ;  e(t)  increasing 

J  n  ^ 


w.(t-t')  ~  ;  e(t)  decreasing 

—1  at 


In  typical  CAD  applications,  the  input  e(t)  is  modelled  by  a 

monotone  waveform.  Hence  the  bounds  given  in  (69)  are  sufficient  for 

most  cases.  However  in  the  rare  case  where  the  input  e(t)  has  both 

positive  ^und  negative  slopes,  we  can  still  bound  the  response  y^(t)  in 

a  similar  way.  We  do  so  by  considering  separately  the  time  intervals  where 
de(t)  de(t) 

—7- —  >  0  and  — -  <  0  when  we  perform  the  convolution  operation  (for 

at  at 

details,  see  Appendix  E  of  [1]). 


b.)  Non-Step  Incut  (Arbitrary  Initial  Equilibrium) 


Suppose  that  we  desire  bounds  on  the  response  y^(t)  to  a  non-step 
input  e (t)  when  we  are  starting  from  some  non-zero  equilibrium. 


Define : 


y(0)  [ay.(t) 


We  can  derive  bounds  for  the  response  Ay^(t)  to  the  non- step  input 
Ae(t)  by  the  convolution  procedure  discussed  in  part  (a),  since  the 
(A  y,  Ae)  system  is  starting  from  a  zero  equilibrium.  Just  add  y^  (0) 
to  these  bounds  in  order  to  obtain  bounds  on  y^(t) : 


^i 

(t)  _<  y^{t)  _<  y^  (t) 

where , 

f 

^^(t)  =  y^(0)  +  Ay^(t) 

* 

yi(t)  »  y^(0)  +  A^^(t) 

(71) . 


The  justification  for  this  is  that  the  ( A y,  A e)  system  is  described 
by  the  same  differential  equation  as  the  (y,  e)  system: 

A‘y(t)  =  y{t) 

=  [RC]  ^  (xe(t)  -  y(t))  ;  see  (9) 

=  [RC]“^  (x(Ae(t)  +  e(0))  -  (A  y(t)  +  y(0))) 

■  [RC]  ^  (x  Ae(t)  +  xe(0)  -  A  y(t)  -  y(0)) 

t _ t 

cancel,  since  y(0)  =  xe(0) 

-  [RC]'^  (X  Ae(t)  -  A  y(t))  (72). 


c. )  Example :  Saturating  Ramp  Input 

The  non-step  input,  e(t),  considered  in  this  example  will  be 


called  a  "saturating  ramp"  since  it  starts  from  zero  volts  at  time  t=0 


and  rises  linearly  until  it  saturates  at  volts  at  time  t  = 

(see  Fig.  17) . 


e  (t) 


Fig.  17:  Saturating  Ramp  Input 

Furthermore,  in  this  example  the  bounds  calculations  will  be  done  under 
the  assumption  that  the  ramp  rises  sharply  compaured  to  the  characteristic 
time  constants  of  the  RC  networ)i  being  driven.  Specifically,  we  shall 
assume: 

0  <  T  ^  <  min  { (T_  -  T_  ) ,  (T  -  T  )}  (73). 

sat  —  D.  R.  P  D 

11  1 

[There  is  nothing  fund^unental  about  the  assumption  in  (73).  Upper  and 

lower  bounds  for  any  choice  of  T  can  be  calculated  by  the  procedure 
— *-  sat 

outlined  in  this  section;  the  only  difference  being  that  the  convolution 
integrals  in  (69)  could  assume  a  different  form.] 

In  the  following  bounds  calculations  for  the  zero-state  response 
at  node  ^  to  the  saturating  ramp  input  e(t),  we  will  ma)ce  use  of  the 


following  edsbreviation: 
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VIII.  Additional  Computational  Burdens  Arising  from  the 
Presence  of  Resistive  Paths  to  Ground 

The  formulas  for  the  time  constants  T  ,  T  ,  and  T  indicate  that 

K .  D ,  P 

1  1 

the  information  required  in  order  to  compute  bounds  for  leadcy  RC  meshes  is 
as  follows: 

1)  The  i  row  (or  column)  of  the  resistance  matrix  R. 

2)  The  diagonal  elements  of  the  resistance  matrix  R. 

3)  All  elements  of  the  vector  x  (i.e.,  the  equilibrium  voltage 
distribution  for  the  entire  network  when  a  1-volt  source  is 
applied  at  the  input) . 

Items  1)  and  2)  above  do  not  represent  much  of  em  additional 
burden  over  non- leaky  RC  meshes  that  happen  to  contain  resistor  loops.  If 
a  non-leaky  RC  mesh  contains  resistor  loops,  it  is  already  difficult  to 
calculate  the  entries  of  the  R  matrix.  The  presence  of  grounded  resistors 
does  not  madte  this  already  difficult  problem  that  much  harder.  If, 
however,  the  RC  network  under  consideration  is  a  leaky  RC  tree,  than  items 
1)  and  2)  post  a  clear-cut  extra  computational  burden.  The  reason  for 
this  is  that  the  entries  of  the  R  matrix  will  require  some  computation: 
we  cemnot  read  them  off  by  inspection  as  we  can  with  a  non-leaky  RC  tree 

(IJ. 

Item  3)  above  represents  a  clear-cut  extra  computational  burden 
regardless  of  what  the  underlying  non-leaky  RC  mesh  looks  like.  For  a 
non-leaky  RC  mesh,  all  equilibrium  voltages  are  equal.  If  resistive 
paths  to  ground  are  present,  then  some  computation  is  required.  However, 
as  the  following  lemma  shows,  the  extra  computation  for  the  equilibrium 


Proof  of  Lemma: 


Assume  the  input  voltage  source  is  a  step  of  magnitude  v^  volts  at 
time  t=0,  and  that  the  network  is  starting  from  a  zero  equilibrium  condition. 
The  differential  equation  for  the  RC  network  is. 


[RC]v(t)  =  V  -  v(t) 
- 


(85) 


Taking  the  k'"^  row  at  time  t=0^,  we  have: 


Ir^Ov  (0^  -  v^(0^ 

D  eq 

+ 


(86) 


But  v,.{0  )  =  0  for  all  k,  and  v^  (0  )  =  0  for  all  j  except  j=l.  So  (86) 
k  3 

becomes , 


(87)  . 


eq 


V 

We  identify  ^  ^  W~  current  flowing  into  node  at  time 


4* 

t=o  , 


\i  ■  r  ■  % 

s  eq 


(88)  , 


V  \l 


V  R 

o  s 


(89)  . 


This  proves  the  lenna. 


This  lemma  shows  that  computing  for  all  k  is  equivalent  to 
computing  Rj^j^  for  all  k  (i.e.,  the  first  column  (or  row)  of  the  resistance 
matrix  R)  However,  the  first  column  (or  row)  of  R  is  exactly  what  would 
be  needed  if  we  wanted  to  bound  the  response  at  node  .  Hence,  the 
calculation  of  the  equilibrium  voltage  distribution  is  no  worse  computationally 


than  considering  one  more  node  as  an  output  node. 


Note  that  in  the  rare  case  when  there  are  M>1  capacitor  nodes 

connected  directly  to  the  source  by  resistors  R  , . . . ,  R  ,  then  a 

3 «  s 

1  M 

Straightforward  extension  to  the  preceding  lemma  holds : 


M  R^. 

I  r 


V  nodes  k 


(90)  . 


This  meams  that  calculation  of  the  equilibrium  voltage  distribution  is 
equivalent  to  con^iuting  the  first  M  columns  (or  rows)  of  R,  which  is 
computationally  equivalent  to  considering  nodes([l)  , . . .  ^^as  additional 
output  nodes  to  be  bounded. 

The  computational  problem  now  consists  of  finding  the  first  row 
(or  column)  of  R,  the  i  row  (or  column)  of  R  for  each  output  node  ^ , 
and  the  diagonal  elements  of  R.  Unfortunately,  the  R  matrix  is  "global" 
in  character  and  is  difficult  to  compute  if  there  are  resistor  loops  or 
resistive  paths  to  ground.  However,  the  conductance  matrix  G  =  R~^ 
is  "local"  in  character,  sparse  (i.e.,  has  memy  zero  entries),  and  can 
be  determined  virtually  by  inspection  even  when  there  are  resistor  loops 
or  grounded  resistors  present.  So  one  approach  is  to  determine  G 
directly  ^u\d  then  invert  it  to  find  R.  This  approach,  however,  raises 
some  doubts  because  it  finds  all  of  R  instead  of  just  the  entries  of  R 
that  are  actually  needed.  An  alternate  approach  is  given  below. 


Alternate  Approach: 

In  this  alternate  approach,  we  only  calculate  the  entries  of  the 
R  matrix  that  are  needed.  If  the  entry  R..  of  R  is  desired,  the  conductance 


matrix  is 


"personalized"  to  nodes  and  by  eliminating  nodes  in  the 
original  resistive  subnetwork  as  much  as  possible  without  destroying  the 
topology  of  node  (T)  with  respect  to  node  ^ .  We  will  illustrate  this 
method  with  an  example  calculation. 


liJ  in  in 


Fig.  19:  Resistive  Subnetwork  for  Example  Calculation  of 
(all  resistors  are  1-ohm) 

The  resistive  subnetwork  for  the  example  calculation  is  shown  in  Fig.  19. 

We  wish  to  calculate  the  entry  R^j^  of  the  R  matrix  that  describes  this 
resistor  subnetwork.  The  first  step  is  to  “personalize"  the  subnetwork 
to  nodes  and  ^  by  eliminating  as  meuiy  nodes  as  possible  (see  Fig.  20) 


in  2n 


Fig,  20:  Resistive  Subnetwork  from  Fig.  19, 
Personalized  to  Nodes  and  ^ 


Note :  We  will  use  a  prime  to  distinguish  the  subnetwork  in  Fig.  20  from 
the  subnetwork  in  Fig.  19. 


The  key  point  to  notice  about  this  node  elimination  procedure  is 

that. 


(91) 


Hence,  we  can  determine  from  the  much  smaller  conductance  matrix  G 

I 

(rather  than  G) ,  either  by  directly  inverting  G  or  by  using  cofactors: 


"ik  ■  "ik  ■  '■32  ■  =.  >I 


(92) 


The  conductance  matrix  G  can  be  read  off  quickly  from  Fig.  20  as; 


G  = 


ii  -1  .1 

6-^3 


-1 


1  i  n 

3  '  2  6 


(93) 


Using  the  cofactor  method,  we  have: 


This  result  of  R.,  =  —  Q  is  correct,  and  can  be  verified  by  going  back  to 
ik  9 

Fig.  19  without  considering  Fig.  20  at  all. 

In  summary,  this  approach  suggests  a  tradeoff.  Instead  of  doing 
one  calculation  with  the  large  original  conductance  matrix  G,  you  can 
do  several  calculations  with  smaller  "personalized"  conductance  matrices  G 
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Appendix  A:  Proof  of  Facts  1  and  2 


The  proof  of  Fact  1  depends  on  additional  structure  of  the 
conductance  matrix  G  alluded  to  in  Section  III. 

Lemma:  Assuming  all  resistors  in  a  given  resistor  subnetwork  R  are 
positive,  the  conductance  matrix  G  which  describes  R  has  the 
following  properties: 

1)  G  >  0  Vn  ;  (i.e.,  all  diagonal  elements  are  positive) 

nn 

2)  G  <0  Vm^^n  ;  (i.e.,  all  off-diagonal  elements  are 

mn  — 

non-positive) 


Proof  of  Lemma: 

Apply  a  1-volt  source  to  node  of  the  resistor  subnetwork  R. 
Ground  the  datum  and  all  external  nodes  other  than  . 

The  Thevenin  equivalent  circuit  looking  into  node  ^  is  a  positive 
resistor  to  ground  •  Since  we  are  applying  a  1-volt  source  to 

node  ^  ,  positive  current  flows  into  node  G^^>0. 

The  Thevenin  equivalent  circuit  looking  into  any  other  node 
(m^n)  is  either  a  positive  resistor  to  ground  occurs  when  node 

is  not  connected  to  node  through  a  single  resistor) ,  or  is  a 

* 

positive  resistor  in  series  with  a  positive  voltage  source 

occurs  when  node  is  connected  to  node  (n)  through  a  single  resistor) . 


It  is  intuitively  clear  that  the  open-circuit  voltage  at  node  ^  is 
non-negative.  This  fact  is  a  special  case  of  the  "voltage  minimax  theorem' 
(3;  pg.  778). 


In  either  case,  since  we  are  applying  ground  to  node  ,  there  can  be 


no  positive  current  flowing  into  node  =>  G  ^  0 


Proof  of  Fact  1: 


(See  Section  VI  for  the  statement  of  Fact  1.) 

The  proof  below  is  essentially  the  same  as  for  non-leaky  RC  meshes 
given  in  [2I : 

If  i=j  or  i=k,  (24)  is  satisfied  trivially.  If  j*k,  then  (24) 

2 

becomes  i  since  R  is  symmetric,  and  is  true  because  R  is 

positive-definite  and  hence  the  2x2  principal  submatrix  of  R  consisting 
of  elements  in  rows  i  or  k  that  also  lie  in  columns  i  or  k  must  have  a 
positive  determin2Lnt.  Now  suppose  i,  j,  and  k  are  distinct,  and  consider 
the  resistive  subnetwork  R  obtained  from  the  original  resistive  sub¬ 
network  R  by  open-circuiting  all  external  terminals  of  R  except  (Ti  ,  (J) 
,  and  the  datum.  Taking  terminals  (p  ,  (j^'  ,  (k)  of  R  as  terminals 
® '  (D '  O  respectively  i 
matrix : 


(A-1)  . 


Now  consider  G  »  R  the  conductance  matrix  of  R.  The  off-diagonal 
elements  of  G  are  non-positive  (by  the  preceding  lemma) ,  so  in  particular 
_<  0.  Taking  the  inverse  of  R  using  cofactors,  we  see  that 

-_i 

G.  ■  =  (R  )i,  -  =  [  (cof  R  )/(det  R)  ]  <0.  Since  R  is  positive-definite 


R,  we 

see 

that  R 

has  the  following 

’  R.  . 

11 

R.  . 
ID 

ik 

R.  . 
Dr 

R.  . 
DD 

'S.i 

\d 

\k 

I 


=>det  R  >  0,  it  follows  that  cof  R  <  0.  But  cof  R.,  *  (-1)  (R.  .  R, 

]k  —  jk  11  k:3 

R^iRij).  This  proves  (24). 


Proof  of  Fact  2: 


(See  Section  VI  for  the  statement  of  Fact  2.) 

Assume  that  the  step  transition  is  "up",  i.e.,  v  >  v  .  We 

-eq  -o 

want  to  prove  that  v^  (t)  2.  0'  nodes  j  in  the  leaky  RC  mesh.  (The 
arg\ament  for  a  "down"  step  is  entirely  analogous.) 

Note  first  that  v^ (0^)  >  0  for  every  capacitor  connected  to  the 

source  through  a  single  resistor,  and  v^(O^)  «  0  for  all  other  capacitors. 

^  N 

Therefore  v(0  )  lies  in  the  "first  orthant"  of  IR  (the  generalization 

2 

to  N-dimensions  of  the  "first  quadrant"  of  IR  ) . 

We  claim  that  v(t),  once  inside  the  first  orthant,  can  never  leave. 
To  see  this,  first  differentiate  (12)  to  obtain  the  following  differential 


equation  for  the  evolution  of  v(t) : 


v(t)  =  -[RC]~^  v(t)  =  -[c‘^G]  v(t) 


(A-2)  . 


The  off-diagonal  elements  of  G  are  non-positive  (by  the  previous  lemma) , 
so  the  off-diagonal  elements  of  - (C  ^G]  are  non-negative,  i.e.,  -[C  ^G] 
is  a  "quasi-monotone  increasing"  operator  [10].  If  v(t)  were  to  exit  the 

first  orthant,  it  would  have  to  do  so  by  passing  through  the  boundary, 

* 

i.e.,  one  of  the  axis  planes.  Thus  there  would  exist  an  instant  t  and 

*  * 

a  capacitor  k  such  that  v  (t  )  »  0  although  v.  (t  )  >  0  Vj  k.  But 

^  D  ~ 

from  (A-2),  this  implies: 


■ 
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)  =  -  I  V  (t  )  >  0  (A-3) 

* 

So  in  fact  v^(t)  does  not  become  negative  at  t=t  .  In  other  words,  the 
vector  field  (A-2)  never  points  out  of  the  first  orthant  when  evaluated 
on  its  boundary.  (For  a  "down"  step,  ^(o"^)  lies  in  the  negative  "first 


orthant"  and  can  never  leave  it.) 


For  a  more  complete  explanation,  and  consideration  of  limiting 
cases,  see  references  [11],  [12],  and  [13]. 


Appendix  B:  Term-by-Term  Ordering  of  the  Time  Constants 


We  have  the  following  formulas  for  the  time  constants, 


'"r.  “  I 

1  11  Jf 


1  1  k 


(B-1) 


It  is  already  known  that  the  inequality  T  <  T  T  holds  for  any 

R  •  ^  D .  ^  P 
1  1 

leaky  RC  mesh  (just  evaluate  the  state  constraints  (46)  at  time  t=0) . 

However,  it  is  not  immediately  obvious  that  this  inequality  also  holds 

term-by-term  for  each  corresponding  term  in  the  summation  (i.e., 

12  1 

FT”  ^ik  ^k  —  ~  ^ik^'k^k  —  \k^k  each  node  k)  .  The  goal  here  is  to 
establish  this  result.  After  simplification,  the  inequality  we  wish  to 
prove  reduces  to: 


(B-2) . 


Proof  of  (B-2) : 

Starting  from  zero  initial  equilibrium,  apply  an  "up"  step  of  1 
volt  at  time  t»0  to  the  input  of  any  leaky  RC  mesh.  Note  that  v  =  0  and 


V 


X.  Evaluating  (26)  at  time  t»0,  we  obtain: 


Appendix  C:  M-Matrix  Systems 


The  purpose  of  this  appendix  is  to  outline  (from  a  more  abstract 
"systems"  point-of-view)  the  underlying  reason  why  the  waveform  bounding 
results  presented  in  (1]  and  [21  carry  over  to  the  le2dcy  RC  meshes  in 
this  report.  Leaky  RC  meshes  are  similar  to  non-leaky  RC  meshes  [2]  and 
non-leaky  RC  trees  [1]  in  one  fundamental  respect:  They  are  all  linear 
dynamical  systems  with  am  "M-matrix"  as  the  system  matrix  [14].  In  the 
literature,  there  are  many  different  equivalent  statements  adjout  M-matrices 
that  would  all  serve  equally  well  as  definitions.  The  statement  best 
suited  to  our  purposes  is  as  follows. 

Definition:  Let  M  be  a  non-singulau:  square  matrix  of  real 
numbers,  amd  let  P  *  m"^.  Then  M  is  said  to 
be  an  "M-matrixf'  if  all  off-diagonal  elements 
of  M  are  non-positive  and  all  elements  of  P 
are  non-negative. 

In  terms  of  the  RC  networks  in  this  report,  the  M-matrix  is 

M  =  c’^G 

and  the  non-negative  inverse  M-matrix  is 

P  »  RC 


(C-1) , 


(C-2)  , 


as  we  can  easily  verify. 

The  conductance  matrix  G  has  non-positive  off-diagonal  elements 


(by  the  lemma  in  Appendix  A) .  Premultiplying  by  the  positive  diagonal 
matrix  C~^  does  not  change  this,  i.e.,  [C  ^G]  has  non-positive  off-diagonal 
elements . 

The  resistance  matrix  R  is  clearly  non-negative,  assuming  that 
the  individual  resistors  are  non-negative.  Postmultiplying  by  the 
positive  diagonal  matrix  C  does  not  chamge  this,  i.e.,  [RC]  is  non-negative 

Hence  [C  ^G]  is  indeed  eui  M-matrix,  as  is  G  itself.  (In  general, 
the  set  of  M-matrices  is  invariant  under  the  group  action  consisting  of 
right  and  left  multiplication  by  positive  diagonal  matrices.) 

The  network  differential  equations  (12)  can  now  be  viewed  more 
abstractly  as  a  linear  dynamical  system  with  a  (negative)  M-matrix  as 
the  system  matrix: 


v(t)  =  -M  v(t)  +  M  V 

,  .  -  .eq 


We  caul  successfully  bound  the  step  response  of  the  abstract  system, 
(C-3) ,  because  both  Fact  1  (appropriately  modified)  and  Fact  2  from 
Section  VI  apply  to  the  abstract  system. 


Fact  1-C:  Let  P  be  a  non-singular  NxN  square  matrix  which 
is  the  inverse  of  an  M-matrix.  Then  P  has  the 
following  structure: 


p .  . p.  .  >  p,  . p .  .  for  auiy  i,j,k  g{1,..., 
k]  —  ki  1] 


(C-4) . 
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Proof  of Fact  1-C: 


If  i=j  or  i=k,  (C-4)  is  satisfied  trivially.  If  j=k,  {C-4)  becomes 

^i^^kk  -  ^ki^ik  -^ki^^ik  general;  the  inverse  M- 

matrix  P  is  not  necessarily  symmetric) ,  which  is  true  because  any 
principal  minor  of  an  inverse  M-matrix  is  known  to  be  positive  [16;  Cor.  1, 
p.  198]  and  hence  the  2x2  principal  submatrix  of  P  consisting  of  the 
elements  in  rows  i  or  k  that  also  lie  in  columns  i  or  k  must  have  a 
positive  determinant. 

Suppose  i,  j,  and  k  are  distinct,  jmd  consider  the  matrix  P  which 
is  the  3x3  principal  submatrix  of  P  consisting  of  the  elements  in  rows  i, 


j,  or  k  that  also  lie  in  columns  i,  j,  or  k.  Since  any  principal  sub¬ 
matrix  of  an  inverse  M-matrix  is  known  to  be  an  inverse  M-matrix  [15,  p.  329], 
*  *_1 

P  IS  an  inverse  M-matrix  and  hence  (P  )^.  <  0.  Taking  the  inverse  of  P 
using  cofactors,  we  have  0  >  =  Kcof  p^j^/(det  P)  ]  ,  and  det  P  >  0 

[16].  Thus  cof  p^j^  =  ^ki^ij  "  ^ii^kj  -  {C-4).  ^ 


Note :  Although  not  proven  here.  Fact  1-C  also  holds  for  a  slightly 
more  general  class  of  matrices  than  inverse  M-matrices.  Going  back  to 
the  more  concrete  domain  of  RC  networks,  if  any  individual  resistor  or 
capacitor  is  equal  to  zero  then  Fact  1-C  still  holds  for  P  =  RC  even 
t.hough  P  is  not  invertible  in  this  case.  However,  P  ^  contained  in  the 
closure  of  the  set  of  inverse  M-matrices.  The  practical  advantage  of 
this  is  that  the  bounding  technique  is  still  valid  for  circuits  with  minor 
specification  errors  (e.g.,  leaving  out  a  capacitor). 
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